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SUMMARY 


Many receptor-mediated endocytic processes are 
mediated by constitutive budding of clathrin-coated 
pits (CCPs) at spatially randomized sites before 
slowly pinching off from the plasma membrane (60- 
100 s). In contrast, clathrin-mediated endocytosis 
(CME) coupled with regulated exocytosis in excitable 
cells occurs at peri-exocytic sites shortly after vesicle 
fusion (~10 s). The molecular mechanism underlying 
this spatiotemporal coupling remains elusive. We 
show that coupled endocytosis makes use of pre- 
formed CCPs, which hop to nascent fusion sites 
nearby following vesicle exocytosis. A dynamic 
cortical microtubular network, anchored at the cell 
surface by the cytoplasmic linker-associated protein 
on microtubules and the LL5B/ELKS complex on the 
plasma membrane, provides the track for CCP hop- 
ping. Local diacylglycerol gradients generated upon 
exocytosis guide the direction of hopping. Overall, 
the CCP-cytoskeleton-lipid interaction demon- 
strated here mediates exocytosis-coupled fast recy- 
cling of both plasma membrane and vesicular pro- 
teins, and it is required for the sustained exocytosis 
during repetitive stimulations. 


INTRODUCTION 


Clathrin-mediated endocytosis (CME) mediates constitutive, 
receptor-mediated endocytosis in all eukaryotic cells. In 
neuron and endocrine cells, rapid release of neurotransmitters 
and hormones is accompanied by massive addition of newly 
fused vesicle membrane to the plasma membrane. Subse- 
quently, coupled endocytosis rapidly initiates to clear the 
plasma membrane for the next round of fusion (Neher, 2010) 
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and maintain the balance of cell surface membrane (Haucke 
et al., 2011). 

CME initiates with the budding of clathrin-coated pits (CCPs) 
that recruit numerous proteins, which bend surface membrane, 
interact with specific cargo proteins, facilitate clathrin coat as- 
sembly, and mediate final membrane scission (McMahon and 
Boucrot, 2011). In non-excitable cells, CCPs constitutively bud 
at spatially random sites before slowly pinching off from the 
plasma membrane with cargo proteins to be retrieved (60- 
100 s) (Ehrlich et al., 2004; Loerke et al., 2009). In neuron and 
endocrine cells, membrane depolarization evokes massive 
exocytosis followed by coupled fast CME (time constants, 
3-10 s) (Artalejo et al., 2002; He et al., 2008; Jockusch et al., 
2005; Wu et al., 2009; Zhu et al., 2009). Moreover, regulated 
CME occurs mostly at the spatially confined peri-active zone 
after stimulation in neurons (Marie et al., 2004; Roos and Kelly, 
1999; Wahl et al., 2013). How the slow, spatially randomized 
CME is synchronized to regulated exocytosis in time and space 
in excitable cells remains enigmatic. 

Using dual-color total internal reflection fluorescence (TIRF) 
microscopy to monitor discrete exocytic and endocytic 
events simultaneously, here we visualize elementary exocy- 
tosis-coupled CME in insulin-secreting INS-1 cells. In addition 
to direct hits of exocytosis onto sites close to pre-formed 
CCPs or de novo formation of CCPs at fusion sites after exocy- 
tosis, we reveal an unexpected hopping of pre-formed CCPs to 
nearby nascent fusion sites. These hopping CCPs move on the 
plasma membrane, with cargo proteins such as Syt VII, to fusion 
sites before endocytosis. Cortical microtubules underneath the 
plasma membrane provide tracks for the hopping of CCPs, 
and transient diacylglycerol (DAG) gradients at fusion sites 
convey the spatial orientation of the movement. Overall, we pre- 
sent a concept that lateral hopping of CCPs with directionality 
precisely couples CME with exocytosis in space and time. In 
conjunction with membrane-cytoskeleton interaction and DAG 
gradient, hopping of CCPs may represent an efficient method 
to actively recycle vesicular components and clear fusion sites 


in excitable cells. 
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The Elementary Coupling between CME and Exocytosis 
We labeled VAMP2 with the fluorophore pHluorin, which remains 
quenched by acidic luminal pH until a vesicle fuses. At a cell foot- 
print, a sudden focal brightening indicates a fusion pore opening; 
and the ensuing circular spread of the pHluorin signal signifies 
VAMP2 proteins diffusing away from the fusion site (Figure S1A; 
Allersma et al., 2004). To simultaneously detect endocytosis, we 
tagged the light chain of clathrin (LCa) with DsRed (Ehrlich et al., 
2004; Merrifield et al., 2002), so that the budding off of a CCP 
could be recorded as the sudden disappearance of an LCa 
punctum (Figures S1A and S1B; Ehrlich et al., 2004; Loerke 
et al., 2009; Merrifield et al., 2005). Stimulation with high glucose 
and high potassium evoked massive exocytosis and the disap- 
pearance of LCa puncta over the entire cell footprint. Visual 
examination of individual VAMP2-pHluorin and LCa image pairs 
revealed that, within an observation window of 150 s, discernible 
LCa puncta transiently existed in close proximity to 38% of 
fusion sites (149 + 7 nm, 151 events). 

To establish an unbiased correlation between exocytosis and 
LCa puncta, we aligned and averaged time-lapse image segments 
of all fusion events and the corresponding LCa signals, using 
the sharply defined pHluorin peak signals to guide the spatio- 
temporal alignment. Ensemble averages revealed that exocytosis 


10 20 30 40 50 
Time (s) 


500 1000 
Length (nm) 


Hopping CCP departure time 


Figure 1. Pre-formed LCa Puncta Hop 


1.04 toward Nascent Fusion Sites for the Exo- 
1.02 3 Endocytosis Coupling 

k (A) Spatiotemporal association of clathrin with 
1.00 exocytosis in INS-1 cells co-transfected with 


VAMP2-pHluorin and LCa-DsRed. (Left) Time- 
lapse images show the average VAMP2-pHluorin 
signal (467 fusion events) at fusion sites (top) and 
the corresponding average LCa signal at fusion 
sites (middle) and at random sites (n = 467, bot- 
tom). (Right) Time courses of normalized fluo- 
rescence signals from average images on the left 
(top) and spatial profiles of the average VAMP2- 
pHluorin and LCa signals at fusion sites (bottom). 
Inset is the overlaid image. 

(B) An example of a hopping CCP (arrowheads, see 
also Movie S1). (Left) Time-lapse images show 
VAMP2-pHluorin (top) and LCa signal (middle) at 
a fusion site, at 3-s intervals. Kymographs (bottom) 
are 150 s and 1.3 um high. (Right) Time courses 
show normalized VAMP2-pHluorin (green) and 
LCa (red) signals at the fusion site. The dotted line 
represents the time course of normalized LCa 
signal at the original site before hopping, and 
departure time (t) represents the latency of a CCP 
departure after a fusion event. 

(C) Normalized histogram of all exocytosis- 
coupled CCP departure times at fusion sites is 
shown (n = 215). 

(D) Normalized histogram of departure times of 
hopping CCPs at fusion sites (n = 74) is shown. 
Scale bars, 500 nm. See also Figure S1 for other 
types of CCP recruitments at fusion sites and 
hopping CCP dynamics analysis and Figure S2 for 
coupled recruitment of Dyn1 at fusion sites. 
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(n = 467) evoked a focal LCa accumulation that overlapped with 
the VAMP2-pHluorin peak in space (Figure 1A). The departure 
time of an LCa punctum from the fusion site (t) was calculated 
as the interval between the moments that the VAMP2-pHluorin 
signal reached a peak and the sudden disappearance of the 
LCa punctum (Figure 1B). We found that the histogram of depar- 
ture time followed a two-Gaussian distribution, with the centers 
at 6.9 and 28.9 s (Figure 1C). Thus, direct visualization reveals 
fast and slow exocytosis-coupled clathrin recruitments. 

In addition to the endocytosis of CCPs, some of the LCa 
puncta disappearance events might have been due to the non- 
productive disassembly of clathrin clusters (Aguet et al., 2013; 
Ehrlich et al., 2004; Loerke et al., 2009). To exclude that possibil- 
ity, we investigated the spatiotemporal correlation of exocytosis 
with labeled dynamin-1 (Dyn1) (Figure S2A), a mechanochemical 
enzyme required for the membrane scission (Aguet et al., 2013; 
Ferguson and De Camilli, 2012). The average Dyn1 signal 
at fusion sites rose immediately upon fusion (Figure S2B). The 
subsequent Dyn1 departure displayed fast and slow Gaussian 
kinetics, with average times of 10.2 and 23.1 s, similar to those 
of LCa puncta (Figure S2C). These results suggest that the 
recruitment of LCa to fusion sites and the ensuing fast and 
slow departures represent coupled CME rather than non-pro- 
ductive events. Collectively, CME is exquisitely orchestrated 
and spatiotemporally coupled with exocytosis. 
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By analyzing individual events from time-lapse movies, we 
identified three elementary types of coupled LCa recruitment. 
Some LCa puncta resided closely to future fusion sites (~31%, 
Figure S1A), as if the fusing vesicles directly hit upon pre-existing 
CCPs, whereas the rest involved LCa puncta either formed de 
novo at fusion sites (~35%, Figure S1B) or that were rapidly 
recruited in toto (hopping) from nearby regions to fusion sites 
(~34%, Figure 1B; Movie S1) after exocytosis. On average, the 
disappearance of pre-existing LCa puncta from the fusion sites 
occurred with a mean latency of 13.1 + 1.6 s (n = 80) after exocy- 
tosis and accelerated coupled CME 5- to 7-fold over classical 
CME (60-100 s in non-excitable cells) (Ehrlich et al., 2004; Loerke 
et al., 2009; Merrifield et al., 2005). This result agrees with the 
concept of a readily retrievable pool, in which pre-assembled 
vesicle membrane proteins wait on the pre-synaptic membrane 
to be rapidly retrieved after exocytosis (Hua et al., 2011). The 
coupled endocytosis of CCPs formed de novo at fusion sites 
was relatively slow (Figure S1B). About 30% of these CCPs did 
not disappear during the observation window of 150 s, and the 
other 70% disappeared from the fusion sites at 35.6 + 3.1 s 
(n = 61) after exocytosis, which was three times slower than 
pre-existing LCa puncta (p < 0.001). Thus, the exocytosis- 
coupled de novo CCP formation and ensuing endocytosis 
operate on a timescale comparable to classical CME (Ehrlich 
et al., 2004; Loerke et al., 2009; Merrifield et al., 2005). 

The discovery of fusion-site-directed, sudden translocation of 
pre-formed LCa puncta was a surprise, because the CCP might 
be thought of as an immobile structure where membrane invag- 
ination, protein recruitment, and vesicle fission occur at the same 
location (Ehrlich et al., 2004; Loerke et al., 2009; McMahon and 
Boucrot, 2011; Merrifield et al., 2005). In search of the possible 
mechanism, we first quantified the movement of LCa puncta. 
On average, an LCa punctum traversed 368 + 32 nm at speeds 
ranging from 36 nm/s to 1.08 um/s (n = 74, Figures S1C and 
S1D), via smooth and sometimes linear trajectories, to reach 
nascent fusion sites (Figure 1B). Once initiated, this hopping 
was completed within ~2.2 + 0.2 s. These characteristics are 
distinct from the long-range (Keyel et al., 2004), dynactin-depen- 
dent fast lateral movement of endosomal clathrin (Zhao et al., 
2007). On average, LCa puncta initiated hopping ~15 s after 
exocytosis, with a broad latency ranging from immediate (2 s, 
25t percentile) to delayed (21 s, 75" percentile, Figure S1E), 
indicating that hopping LCa puncta participate in both the fast 
and the slow coupled CME. The average departure time was 
28.2 + 2.2 s (exclusive of events outlasting the observation win- 
dow, Figure 1D), which was ~15 s slower than that of pre-exist- 
ing CCP-mediated CME but ~7 s faster than that of CCPs 
formed de novo at fusion sites. Combined, the direct targeting 
of exocytosis onto sites close to pre-existing CCPs and the 
lateral hopping of LCa puncta toward fusion sites provide a 
long-sought local mechanism of rapidly supplying endocytic 
structures during exocytosis. 


Pre-formed CCPs Hop with Cargo Proteins on the 
Plasma Membrane toward Nascent Fusion Sites for 
Protein Retrieval 

It has been postulated that clathrin can be a component of the 
fusing granules in retinal bipolar neurons (Pelassa et al., 2014), 
which may contribute to the hopping LCa puncta and couples 


CME with exocytosis. We used a fast incidence angle-scanning 
and azimuthal-averaging TIRF microscope to perform super- 
resolution reconstruction in the z axis to directly test that possi- 
bility (Boulanger et al., 2014; Zong et al., 2014). As shown in 
Figure 2A, the sudden approaching of a secretory vesicle 
(labeled with VAMP2-EGFP) along the z axis was not colocalized 
with any LCa puncta before collapsing into the plasma mem- 
brane. After the completion of the exocytosis, a nearby pre- 
formed LCa punctum, seen only at the plane 0-50 nm beneath 
the plasma membrane, moved toward the fusion site and 
docked there. After a delay of ~39 s, the shallow clathrin struc- 
ture started to invaginate and appeared in the plane that was 
100-150 nm beneath the plasma membrane. Finally, it sequen- 
tially disappeared from both planes, indicative of internalization 
of a clathrin-coated vesicle (Figures 2A-2C). This result indicates 
that hopping LCa puncta represent CCPs associated with the 
plasma membrane. 

Next we explored whether hopping CCPs contain cargo 
proteins and other membrane components. Synaptotagmin VII 
(Syt VII) is a vesicular Ca2* sensor mediating glucose-induced 
insulin secretion (Gustavsson et al., 2008), and it needs to be 
retrieved from the plasma membrane by CME after vesicle 
exocytosis (Haucke and De Camilli, 1999). By photo-activating 
Syt VII-PA-GFP/PA-mCherry under TIRF illumination, we visual- 
ized all Syt VII proteins on the plasma membrane as well as a 
small fraction of vesicular Syt VII close to the plasma membrane; 
furthermore, by briefly incubating cells with a cell-impermeable 
dye (30 s), we visualized a fraction of CLIP-Syt VII proteins local- 
ized specifically on the plasma membrane. With both protocols, 
we observed discrete patches of Syt VII left on the cell surface 
after stimulation (Figure S3A), similar to previous observations 
in synapses (Hua et al., 2011; Willig et al., 2006). Most Syt VII 
patches co-localized and laterally moved with LCa puncta on 
the plasma membrane before ultimately disappearing altogether 
(Figures 2D, S3B, and S3C; Movie S2). In addition, hopping Syt 
VII puncta traveled to vesicle fusion sites (Figure 2E) and disap- 
peared there 4.8 or 17.2 s after exocytosis (Figure 2F), qualita- 
tively matching the kinetics of LCa puncta. This demonstrates 
that clathrin hops to fusion sites along with Syt VII to be internal- 
ized together. We also proved that these LCa puncta slid with 
membrane curvatures. FCHO2 binds and stabilizes shallow cur- 
vatures on the surface membrane (Shen et al., 2012). Most LCa 
puncta co-localized with FCHO2 (71.4% + 0.9%, seven cells), 
and many LCa puncta laterally shifted on the plasma membrane 
with FCHO2 (Figure 2G; Movie S3). The laterally shifted FCHO2 
puncta also moved to the fusion sites defined by the transient 
brightening of the large-dense-core vesicle marker NPY- 
pHluorin (Figure S4A). Finally, AP2, an adaptor protein critical 
for the cargo sorting of CCPs, also was found to laterally move 
on the plasma membrane toward the fusion site after vesicle 
exocytosis (Figure S4B). Taken together, we conclude that mem- 
brane-protein assemblies in toto, i.e., CCPs, hop to the fusion 
sites to mediate recycling of vesicular proteins. 


Cortical Microtubules Provide Tracks for the Hopping 

of CCPs 

The cortical cytoskeleton, which regulates diffusion of mem- 
brane receptors on the plasma membrane (Jagqaman et al., 
2011), interacts with CCPs (Montagnac et al., 2013). Thus, we 
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Figure 2. Pre-formed CCPs Hop with Cargo Proteins on the Plasma Membrane toward Nascent Fusion Sites to Be Internalized 

(A) A representative example of a vesicle fusion event and correlated hopping LCa puncta event. INS-1 cells were transfected with VAMP2-EGFP and LCa- 
DsRed, stimulated with 70 mM KCI and 15 mM glucose. (Left) The cartoon shows different clathrin-coated structures visualized by the super-axial resolution 
incidence angle-scanning and azimuthal-averaging TIRF microscope. (Right) Time-lapse images show the VAMP2-EGFP (green, top) and the LCa signal (red, 
middle, LCa[s]) at the plane 0-50 nm beneath plasma membrane and the LCa signal at the plane 100-150 nm beneath plasma membrane (blue, bottom, LCa|[d)). 
Montages are shown at 32-s intervals. Arrows indicate the moving CCP. 

(B) Time courses of normalized VAMP2-EGFP (green) and LCa (red) signals at the fusion site from the plane 0-50 nm beneath plasma membrane. The dotted line 
represents the time course of normalized LCa signal at the original site before hopping. 

(C) Time courses of normalized LCa signals at the fusion site from the planes 0-50 nm (red) or 100-150 nm (blue dotted) beneath plasma membrane are shown. 
(D) A representative example of an LCa/Syt VII co-cluster (arrowheads, top) that migrated on the plasma membrane before disappearing altogether (see also 
Movie S2). CLIP-Syt VII was labeled with non-permeable CLIP-Surface 488. Montage images are at 19.8, 23.1, 26.4, 33.0, 36.3, 42.9, 46.2, and 49.5 s. 
Kymographs (bottom) are 3.2 um high. Red arrow indicates the time point of linear movement; arrowheads indicate sites where intensities are plotted on the right 
graph. Time courses of normalized Syt VII (green) and LCa (red) signals of the event are shown (right). Inset describes the labeling strategy. 

(E) Arepresentative example of a Syt VII cluster that moved to and then disappeared at the fusion site after exocytosis. Time-lapse images of VAMP2-pHluorin 
(top) and of Syt VII-PA-mCherry signal (middle) at a fusion site are shown at 1.5-s intervals. Kymographs (bottom) are 2.4 um high. Green arrow, time point of 
vesicle fusion; red arrow, the starting time point of Syt VII movement. Time courses of normalized fluorescence signals are shown (right). Inset describes the 
labeling strategy. 

(F) Histogram of Syt VII cluster departure times (t) at fusion sites is shown (n = 144). 

(G) Arepresentative example of a hopping CCP (LCa-EGFP, arrowheads) that co-localized with an FCHO2-mCherry punctum, migrated on the plasma membrane 
before disappearing altogether (see also Movie S3). Montage images (top) are at 3-s intervals. Kymographs (bottom) are 1.7 um high. Red arrow marks the 
initiation of LCa/FCHO2 movement. Scale bars. 500 nm (A) and 1 um (D, E, and G). See also Figure S3 for the co-localization and co-movement of plasma 
membrane-associated Syt VII with CCPs and Figure S4 for the hopping of other adaptor proteins of CCP toward nascent fusion sites labeled by different markers. 


hypothesized that the cytoskeletal network immediately with nocodazole, but not by depolymerization of filamentous 


beneath the plasma membrane might provide tracks for the 
CCP hopping. By labeling microtubules with EGFP-Tubulin, 
we observed that a CCP moved along a cortical microtubule 
filament and stalled at the intersection of two filaments; as 
this microtubule re-oriented, the CCP went along with the 
sliding filament (Figure 3A; Movie S4). To quantitatively correlate 
CCP motions with cytoskeleton network, we tracked trajec- 
tories of all CCPs and identified ~10% linearly moving CCPs, 
using methods proposed previously (Jagaman et al., 2008, 
2011; Figures 3B [inset] and 3C). These linear CCP movements 
were preferentially overlapped with cortical microtubules more 
than with actin filaments (labeled with Lifeact-EGFP, Figure 3B). 
Consistently, the proportion of CCPs that exhibited linear 
motion was reduced by disruption of the microtubular network 
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actin with either cytochalasin D (CytoD) or latrunculin A (Figures 
3C and S5A). This is different from the actin-dependent lateral 
movement of clathrin structures reported previously (Gaidarov 
et al., 1999). We propose that CCP hopping is a subset of linear 
CCP movements, since nocodazole treatment reduced both the 
overall CCP recruitment and the hopping of CCPs to fusion 
sites, whereas CytoD treatment reduced neither (Figures 3D 
and 3E). In addition, inhibition of kinesin 5, dynactin, or myosin 
5 all failed to reduce the linearly motile population of CCPs 
(Figure 3C). Given the relatively slow speed of CCP hopping 
(Figure S1D), these results suggest that linear CCP motion is un- 
likely to be regulated by classical motor proteins, but may result 
from an association of CCPs with dynamic cortical microtubules 
(Montagnac et al., 2013). 
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Figure 3. CLASP Mediates the Anchoring of Cortical Microtubules on the Surface Membrane, which Regulates Both Lateral Movement of 
CCPs and Regulated Exocytosis 

(A) Time-lapse images show a mobile CCP (LCa-DsRed, arrowheads) moving along a dynamic microtubule filament (EGFP-Tubulin, white line) and then being 
redirected by the sliding microtubule. See also Movie S4. 

(B) The ratios of linear CCP trajectories overlapped with actin filaments (Lifeact-EGFP, 3,817 trajectories from six cells) or microtubules (EGFP-Tubulin, 8,262 
trajectories from nine cells, p < 0.001). Inset shows representative linear CCP trajectory (blue) overlaid on EGFP-Tubulin (white). 

(C) Percentages of linearly moving CCPs in LCa-DsRed-expressing control cells (n = 14); in LCa-DsRed-expressing cells incubated with CytoD (20 uM, 20 min, 
n = 16), latrunculin A (10 uM, 30 min, n = 27), or nocodazole (10 uM, 2 hr, n= 7, p < 0.001); or in LCa-DsRed-expressing cells co-transfected with KHC-DN (n = 4), 
p50 (n= 7, p < 0.01), or myo5a-DN (n = 9, p < 0.05), which inhibit kinesin, dynein, or myosin5, respectively. Unless specified otherwise, all statistical differences are 
between the control cells and the designated treated cells. 

(D) The percentage of exocytosis associated with CCP recruitment in control cells (n = 6) and cells pre-treated with CytoD (n = 5) or nocodazole (n = 5, p < 0.001) is 
shown. 

(E) The percentage of hopping CCP among the total CCP recruitment at fusion sites in control cells (n = 6) and cells pre-treated with CytoD (n = 5) or nocodazole 
(n = 5, p < 0.01) is shown. 

(F) An example of a microtubule filament (TagRFP-T-Tubulin, arrows) sliding back and forth between two pairs of ELKS puncta (GFP-ELKS, arrowheads) upon 
stimulation is shown. See also Movie S5. Simplified cartoons are shown (bottom). 

(G) Microtubule network and their skeletal extraction pictures of microtubule filaments under TIRF illumination in control cells (INS-1 cells transfected with EGFP- 
Tubulin, left) and cells co-transfected with EGFP-Tubulin and mTagBFP2-CLASP2-AC (right) are shown. 

(H) The percentage of linearly moving CCPs in control cells (n = 8) and cells overexpressing CLASP2-AC (n = 9, p < 0.001) is shown. 

(I) The percentage of exocytosis associated with CCP recruitment in control cells (n = 19) and cells overexpressing CLASP2-AC (n = 8, p < 0.01) is shown. 

(J) The membrane capacitance increment evoked by trains of depolarization in control cells (n = 10) and cells overexpressing CLASP2-AC (n = 7, p < 0.05) is 
shown. Scale bars, 1 um (A), 500 nm (B), 2 um (F), and 5 um (G). See also Figure S5A for the disruption of actin filament by CytoD or latrunculin A treatment, Figures 
S5B-S5D for the stimulation-evoked accumulation of LL5B/ELKS clusters on the plasma membrane, and Figure S5E for the capacitance traces in control cells 
and in cells transfected with CLASP2-AC mutant. 


Next we tried to search for proteins anchoring cortical micro- HeLa cells (Lansbergen et al., 2006). In INS-1 cells co-trans- 


tubules to the plasma membrane. Microtubule filaments are 
found in close proximity to the active zone in synapses (Westrum 
and Gray, 1986), and they are attached to the cell cortex via the 
interaction of cytoplasmic linker-associated protein (CLASP) on 
the microtubules with LL5B/ELKS on the plasma membrane in 


fected with TagRFP-T-Tubulin and EGFP-ELKS, we found that 
most cortical microtubules traversed ELKS puncta; after stimu- 
lation, more LL5B/ELKS puncta appeared on the plasma mem- 
brane (Figures S5B-S5D). Correspondingly, segmented fila- 
ments frequently slid between different pairs of ELKS puncta 
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(Figure 3F; Movie S5). In cells transfected with a CLASP2 trunca- 
tion mutant (CLASP2-AC) that is unable to interact with the LL5p/ 
ELKS complex (Lansbergen et al., 2006), cortical microtubules 
were reduced (Figure 3G), and both the linear CCP movement 
and the exocytosis-coupled CCP recruitment at fusion sites 
were inhibited (Figures 3H and 3I). Depolarization-evoked exocy- 
tosis also was inhibited in these cells (Figures 3J and S5E), sup- 
porting a role of CLASP and ELKS in regulated exocytosis as 
shown previously in synapses and endocrine cells (Beffert 
et al., 2012; Haucke et al., 2011; Ohara-Imaizumi et al., 2005). 
Taken together, the dynamic, CLASP-dependent anchoring of 
cortical microtubules to the plasma membrane constitutes an 
exo-endocytosis coupling factor. 


Local, Transient DAG Gradients Convey the Spatial 
Orientation for the Hopping of CCPs 

To initiate and direct CCP hopping along the microtubular track, 
the triggering signal must be exocytosis specific, local, and tran- 
sient. In this regard, depolarization-triggered Ca?* influx that is 
required for exocytosis may activate phospholipase C, which 
hydrolyzes phosphatidylinositol 4,5-bisphosphate to generate 
DAG (Kadamur and Ross, 2013). DAG is a non-bilayer lipid 
with a negative curvature (Carrasco and Mérida, 2007), and its 
accumulation may reduce the elastic energy of membrane defor- 
mation and facilitate lateral CCP movement. Therefore, we used 
the EGFP-tagged tandem C1 domains of protein kinase D (PKD) 
as a DAG indicator (PKD-C1ab-EGFP) (Kim et al., 2011), and we 
tested whether the local DAG gradient could be a triggering 
signal. Global increases in DAG concentration were detected 
in cells stimulated with carbachol (Figure S6). More importantly, 
by spatiotemporal averaging, we found that the local DAG con- 
centration began to increase at the fusion sites prior to exocy- 
tosis, reached a plateau when the fusion pores fully opened, 
and stayed elevated for some time after exocytosis (Figures 4A 
and 4B). Incubating cells with a DAG analog, 1-oleoyl-2-acetyl- 
sn-glycerol (OAG), reduced the surface membrane tension as 
measured by optical tweezers (Figures 4C-4E) and increased 
the speed of linear CCP movement (Figure 4F). In contrast, inhi- 
bition of phospholipase C with U73122, which blocked DAG gen- 
eration, reduced the speed of linear CCP movement (Figure 4G). 
The speed of linearly moving CCPs was consistently faster close 
to the fusion sites than elsewhere (Figure 4H). Hence, our data 
suggest that the exocytosis-associated local DAG gradient initi- 
ates and directs the hopping of CCPs toward nascent fusion 
sites by reducing the energy barrier against intra-membrane 
movement of CCPs. This also highlights a general, unappreci- 
ated role of membrane microdomains with transient curvatures: 
they may attract curved sub-cellular structures in addition to the 
recruitment of cytosolic proteins with curvatures as proposed 
recently (Galic et al., 2012). 


Exocytosis-Coupled CME Contributes to the Fast 
Recycling of Both Plasma Membrane and Vesicular 
Proteins 

In pancreatic B cells, membrane depolarization triggers CME and 
clathrin-independent endocytosis; both processes are dynamin 
dependent (He et al., 2008). Thus, we used pitstop (von Kleist 
et al., 2011) to block CME and dynasore (Macia et al., 2006) to 
block CME plus clathrin-independent endocytosis. Using repet- 
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itive trains of depolarization in control cells and in cells pretreated 
with either pitstop or dynasore, we examined evoked exocytosis 
and endocytosis. Inhibition of cellular dynamin function with dy- 
nasore eliminated most of the exocytosis-coupled endocytosis 
(Figures 5A, 5B, and 5D). Pitstop pretreatment increased the 
time constant of the fast capacitance decay (control: 0.35 + 
0.04 s, n = 8; pitstop treated: 0.66 + 0.08 s, n= 8, p < 0.01) (Fig- 
ure 5A), and it reduced the membrane retrieval within 4 s after a 
train of depolarization by half (Figure 5D). Combined, CME con- 
tributes to ~50% of the fast membrane retrieval after exocytosis. 
On the other hand, inhibition of CME failed to affect exocytosis 
evoked by the first train of depolarization, but significantly in- 
hibited subsequent exocytosis evoked by the fourth and fifth 
trains of depolarization (Figures 5A-5C). Inhibition of coupled 
endocytosis also reduces the fast release of synaptic vesicles 
evoked by repetitive stimulations delivered at intervals less 
than ~10 s in Calyx synapses (Hosoi et al., 2009). Therefore, 
fast coupled CME plays an important role in sustaining multiple 
rounds of exocytosis during prolonged stimulations in neurons 
and endocrine cells, as proposed previously (Neher, 2010). 

Next we explored how CME contributed to the vesicular 
protein recycling process by studying endocytosis of surface 
VAMP2 protein tagged with fluorogen-activating protein (FAP) 
(Szent-Gyorgyi et al., 2008). Because the fluorogen dye alone re- 
mains dark until binding to the FAP, its fluorescence intensity is 
directly correlated with the amount of VAMP2-FAP protein 
labeled. Using a cell non-permeable «RED fluorogen, we found 
that surface VAMP2-FAP was internalized and transported to 
puncta and Golgi-like structures in the cytosol during 10 min of 
stimulation. In contrast, most recycling of VAMP2 protein after 
stimulation was blocked in cells pretreated with pitstop, while 
dynasore pretreatment only slightly enhanced the accumulation 
of VAMP2-FAP on the plasma membrane than pitstop pretreat- 
ment (Figures 5E and 5F). Therefore, although CME only acts as 
a complementary route of plasma membrane retrieval, it plays a 
dominant role in the recycling of VAMP2 protein on the plasma 
membrane. 


DISCUSSION 


Overall, we have unraveled an exquisite, DAG-guided and micro- 
tubule-assisted mechanism for the spatiotemporal regulation of 
exo-endocytosis coupling in excitable cells. We envision that 
pre-formed CCPs are motile structures on the plasma mem- 
brane. Excitable cells make use of these pre-formed CCPs to 
accelerate exocytosis-coupled CME, both by directly targeting 
them close to fusion sites and by moving them from nearby re- 
gions. Upon stimulation, pre-formed CCPs laterally hop along 
dynamic cortical microtubules mediated by CLASP on the mi- 
crotubules and LL5B/ELKS on the surface membrane, and they 
are driven by a transient, local DAG gradient toward nascent 
fusion sites (Figure 6). This enables CCPs that initially assembled 
at spatially randomized sites to retrieve vesicle membrane pro- 
teins at preferred nascent fusion sites in a timely fashion. 

Now we have shown that exocytosis-coupled CME and cla- 
thrin-independent endocytosis each contributed to ~50% of 
the fast membrane retrieval after depolarization (Figure 5D). In 
contrast, almost all of the VAMP2 internalization after stimulation 
was mediated by coupled CME (Figures 5E and 5F). These data 
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Figure 4. Local DAG Gradient Generated 
upon Exocytosis Directs the Hopping of 
CCPs toward Nascent Fusion Sites 

(A) Local DAG indicated by PKD-C1ab-GFP tran- 
siently accumulated at vesicle fusion sites indi- 
cated by VAMP2-mOrange2. Time-lapse images 
show the average VAMP2-mOrange2 signal (of 
193 fusion events, top) and the corresponding 
average DAG signal (bottom) at fusion sites. Scale 
bars, 500 nm. Time courses of average fluores- 
cence signals are shown (right). 

(B) Surface plots of the mean DAG gradient (color- 
map coded) around fusion sites at times before 
(28 s), during (32 s), and after exocytosis (36 s) are 
shown. 

(C) Arepresentative example of amembrane tether 
formed as the bead was retracted away from the 
cell is shown. 

(D) Examples show tether force versus distance 
from a control INS-1 cell (left) and an OAG-treated 
cell (100 uM, right). Red, the force required for re- 
tracting the bead away from the cell at a speed of 
1 um/s; gray, the baseline force required for mov- 
ing the bead toward the cell at 1 um/s. The differ- 
ence between these two forces is the calibrated 
tether force. 

(E) Surface membrane tension measured by opti- 
cal tweezers in control (n = 32) and OAG-treated 
(100 uM, 15-20 min, n = 32, p < 0.001) cells is 
shown. 

(F) The velocity of linear CCP movement in OAG- 
treated cells as compared to the control. INS-1 
cells were stimulated with high KCI and high 
glucose in absence or presence of OAG (0 uM 
as control, mean velocity of linear CCP move- 
ment, 598 + 14 nm/s, 2,704 trajectories from 
nine cells; 10 uM, 2,724 trajectories from 
nine cells; 100 uM, 4,451 trajectories from twelve 
cells, p < 0.001). 

(G) The velocity of linear CCP movement in cells 
pre-incubated with the phospholipase C inhibitor 
U73122 (10 uM) for 15 min (3,447 trajectories from 
eight cells) as compared to control cells (mean 
velocity of linear CCP movement, 532 + 40 nm/s, 
1,142 trajectories from four cells, p < 0.01) is 
shown. 

(H) Linear movement of CCPs at fusion sites 
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(27 trajectories) compared to that at random, non-fusion sites (control, mean velocity of linear CCP movement, 543 + 6 nm/s, 1,256 trajectories, p < 0.05). See 
also Figure S6 for the global increase in cellular DAG concentration detected by PKD-C1ab-GFP upon carbachol stimulation. 


highlight a selective recycling of vesicular proteins via fast exocy- 
tosis-coupled CME. Conceptually, lateral hopping of CCPs rep- 
resents an active paradigm in retrieving vesicular proteins on 
the plasma membrane, in contrast to the conventional view of 
CCPs as spatially stationary pits that passively collect cargos 
passing by. By moving toward nascent vesicle fusion sites en- 
riched with cargo proteins, shallow CCPs may act as an active 
sorting machinery to recruit vesicular proteins, especially for 
those non-diffusing ones such as phogrin (Tsuboi et al., 2004). 
The merge of shallow CCPs with cargos reduces abortive CCP 
events and facilitates CCP maturation and effective endocytosis, 
as others have suggested (Swan, 2013). Although pre-existing 
CCPs and de novo formation of CCPs at fusion sites share 
some of these advantages, only hopping CCPs can rapidly and 
continuously supply partially assembled endocytic machinery 


to the vesicle release sites. As compared to the control, exocy- 
tosis evoked by repetitive stimulations was gradually reduced 
in cells with compromised CME, and the reduction started to 
reach statistical significance at the fourth train of depolarization 
(Figures 5A-5C). Given that the reduced exocytosis is mostly 
likely to be due to the disruption of fusion sites by prior exocy- 
tosis, our data further support that hopping CCPs may be impor- 
tant for the local, fast, and efficient clearance of vesicular compo- 
nents left at the fusion sites during repetitive stimulations. 

The microtubule-dependent linear movement of CCP is 
intriguing. It may be a passive process through which the sliding 
microtubules push CCPs attached to them toward fusion sites. 
The -tubulin acetyltransferase («TAT1) in the CCPs, which acet- 
ylates and stabilizes the microtubules (Montagnac et al., 2013), 
could be the physical linker between microtubules and CCPs. 
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Figure 5. Inhibition of CME Significantly Reduced Recycling of Plasma Membrane, Vesicular Proteins, and Local Clearance of Fusion Sites 
upon Repetitive Stimulations 

Five consecutive rounds of depolarization stimulations were delivered to control INS-1 cells (n = 10, black) and INS-1 cells pretreated with either pitstop (24 uM, 
10 min, n = 8, blue) or dynasore (80 uM, 10 min, n = 14, green). The time interval between different rounds of depolarization was 5 s. 

(A and B) The average increases and decreases in membrane capacitance evoked by the first depolarization stimulation (A) and by the fourth depolarization 
stimulation (B) are shown. 

(C) Evoked exocytosis from control INS-1 cells and INS-1 cells pretreated with either pitstop or dynasore during five consecutive rounds of depolarization 
stimulations are shown (statistical significances as indicated on the graph). 

(D) Within 4 s after the very first depolarization, evoked endocytosis in control INS-1 cells and INS-1 cells pretreated with either pitstop or dynasore were 
calculated as decreases in membrane capacitance from 6.2 to 10.2 s in (A). 

(E) Representative images show INS-1 cells transfected with VAMP2-FAP (red) treated with 70 mM KCI, 15 mM glucose solution containing DMSO vehicle (0.8%), 
pitstop (25 uM), or dynasore (80 uM) for 10 min. Scale bar, 5 um. 

(F) The recycling efficiency of VAMP2, quantified as a ratio of average fluorescence intensities of the whole cell boundary to that of the cell cytosol, in a cell treated 
with 0.8% DMSO (n = 120), pitstop (n = 144, p < 0.001), or dynasore (n = 132, p < 0.001). The difference in the VAMP2 recycling is also significant between pitstop- 


treated cells and dynasore-treated cells (p < 0.01). 


Alternatively, cortical microtubule filaments could serve as cages 
that convert random CCP mobility into linear motion, similar to 
the linear diffusion of CD36 receptors on the plasma membrane 
due to the cytoskeletal confinement (Jaqaman et al., 2011). 
Further work is needed to delineate these possibilities. Neverthe- 
less, we did observe an increase in the number of ELKS clusters 
on the plasma membrane upon stimulation (Figures SSB-S5D), 
which was accomplished with dynamic reorganization of cortical 
microtubule networks (Figure 3F). Moreover, CLASP-dependent 
anchoring of cortical microtubules to the plasma membrane 
were required for vesicle fusions, linear CCP motion, and CCP 
recruitment at fusion sites (Figures 3H-3J). Thus, the three- 
dimensional cortical microtubules anchored around ELKS- 
defined fusion sites may represent a spatiotemporal coupling 
framework, which increases the probability of cargo proteins 
merging with mobile CCPs and enables efficient clearance of 
plasma membrane at fusion sites. 

Although this study was conducted in an endocrine cell, func- 
tions of the proteins involved (such as ELKS and CLASP) are 
conserved in synapses, suggesting that a conserved mechanism 
also may account for the fast CME critical for the sustained neuro- 
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transmission in synapses (Kawasaki et al., 2011). Moreover, the 
transient increase in DAG on the plasma membrane close to exo- 
cytic sites observed here also may facilitate generation of clathrin- 
independent invaginations observed shortly after depolarization 
in synapses (Kononenko et al., 2014; Watanabe et al., 2013). 
The microtubule-dependent movement of CCPs facilitates active 
sorting of vesicle membrane proteins on the plasma membrane 
toward designated sites to be internalized, which also may partic- 
ipate in other extracellular signal-regulated fine-tuning of cell 
membrane components in other eukaryotic cells. 


EXPERIMENTAL PROCEDURES 


Live-Cell Imaging and Fluorescence Microscopy 

All experiments were performed at 37°C. INS-1 cells were planted on glass 
coverslips and incubated with bath solution containing the following (in mM): 
136 NaCl, 4.2 KCI, 2.4 CaCl, 1.2 KH2PO;4, 1.2 MgSO,;,, 4 glucose, 10 HEPES, 
and 1 L-glutamine (pH 7.4). Individual coverslips were then placed in a metal 
chamber mounted on a heated stage. Unless specified otherwise, INS-1 cells 
were stimulated with a solution containing the following (in mM): 70 NaCl, 70 
KCl, 2.4 CaCl», 1.2 KH2PO;, 1.2 MgSOuz, 15 glucose, 10 HEPES, and 1 L-gluta- 
mine (pH 7.4) to trigger exocytosis and endocytosis. 
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Figure 6. Model of How Fast CME Is Coupled with Regulated 
Exocytosis 

The local DAG gradient steers the hopping of a mobile CCP toward the fusion 
site along tracks provided by cortical microtubule filaments, which also help 
the delivery of secretory vesicles to fuse with the plasma membrane. 


For Figures 1, 2D-2F, and S1-83, a TIRF microscope was equipped with a 
Sensicam EM682KX camera (PCO AG) and a 100x 1.49 numerical aperture 
(NA) oil-immersion TIRF objective (Olympus). Most images were acquired at 
2 Hz, except for Figures 2D, 2E, and S3 (3.3 Hz). Each pixel from the camera 
was 67 X 67 nm. 

A different TIRF setup was used for Figures 2G, 3, 4, 5, and S4—S6. We used 
an Olympus IX81 inverted microscope equipped with a 150x 1.45 NA oil-im- 
mersion TIRF objective (Olympus). Fluorescence signals were acquired using 
an electron-multiplying charge-coupled device (EMCCD, Andor iXon3 888) at 
a sampling rate of 6.7 Hz controlled by MetaMorph (Molecular Devices). Each 
pixel from the Andor camera was 87 x 87 nm. A 473-nm laser was used to 
excite GFP/pHluorin, a 561-nm laser to excite mCherry/RFP/DsRed/ 
TagRFP-T, and a 405-nm laser to excite mTagBFP2. For two-color imaging, 
100-nm fluorescent beads (TetraSpecK Microspheres, Life Technologies) 
were deposited onto a coverslip and imaged as markers for later alignment. 
To visualize surface VAMP2-FAP internalization (Figure 5E), the TIRF illumina- 
tion was adjusted to wide field. 

The fast incidence angle-scanning and azimuthal-averaging TIRF micro- 
scope was built on a commercial microscopy (1X81, Olympus) equipped with 
an apochromatic objective lens (APO 100x/1.65 NA, Olympus) (Figures 2A- 
2C). Two excitation lasers (Sapphire 488LP and Sapphire 561LP, Coherent) 
were collimated and combined by an AOTF and projected onto a two-axis 
MEMS (A1B2.2, Mirrorcle Tech), which enabled fast incidence angle scanning 
or the scanning of a circle at the back focal plane of the objective for azimuthal 
averaging at 1 kHz. The emitted fluorescence was captured by an EMCCD 
camera (iXon 897U, Andor). The whole system was controlled by Labview 
(Zong et al., 2014). Cells were cultured on customized N-LAF21 coverslip 
with an index of refraction of 1.78. To achieve 50-nm axial resolution, sequen- 
tial dual-color TIRF imaging at ten different incidence angles (61.34, 51.45, 
51.59, 51.795, 52.085, 52.54, 53.29, 54.67, 57.72, and 67.38) were recorded 
at a frequency of 0.8 s/vol. Image reconstruction was conducted as reported 
previously (Boulanger et al., 2014). 


Visualization of Surface VAMP2-FAP Internalization 

In experiments aimed to visualize VAMP2-FAP internalization, INS-1 cells were 
transfected with VAMP2-FAP vector and labeled with a bath solution contain- 
ing 70 mM KCI, 15 mM glucose, and non-permeable aRED fluorogen (100 nM, 
SpectraGenetics) for 10 min at 37°C. These cells were transferred to an imag- 
ing chamber without washing and imaged under wide-field illumination. 


Selectively Highlighting or Labeling Surface Syt VII 
Two methods were used to label the surface pool of Syt VII. In experiments using 
Syt VII-PA-GFP/PA-mCherry (Figures 2E, 2F, and S30), cells were first illumi- 


nated with a 405-nm laser for 15 s to photo-activate the Syt VII population close 
to the cell surface and then imaged during stimulation. In experiments using 
CLIP-Syt VII, cells were first incubated in a stimulation solution containing KCl, 
glucose, and non-permeable CLIP-Surface 488 (New England Biolabs) (Figures 
2D, S3A, and 3B; 20 mM) for 30 s, which was used to label the surface pool of Syt 
VII that was transiently inserted into the plasma membrane due to vesicle exocy- 
tosis. Cells were then washed three times in normal bath solution for 2 min before 
imaging under stimulation with a solution containing KCI and glucose. 


Electrophysiology 

Electrophysiological experiments were conducted by standard whole-cell 
recording using an EPC-10 patch-clamp amplifier (HEKA) as previously 
described (Liang et al., 2011). The capacitance was measured using the 
sine + DC mode of the lock-in extension in the Pulse software, with an 800- 
Hz, 20-mV peak-to-peak sinusoid voltage applied onto a DC holding potential 
of —70 mV. In Figure 5, the intracellular pipette solution contained the following 
(in mM): 125 Cs-glutamate, 10 NaCl, 0.05 EGTA, 1 MgClo, and 5 HEPES, with 
pH adjusted to 7.35 using CsOH. Five trains of stimulation were sequentially 
delivered to the control cells or cells pretreated with inhibitors. Each train of 
stimulation (Cm5+8) consisted of five 50-ms and eight 500-ms depolarizations 
(from —70 to 0 mV) to evoke release from the immediately releasable pool (IRP) 
and the readily releasable pool (RRP) of vesicles, respectively. The inter-pulse 
interval was 100 ms (Merrins and Stuenkel, 2008) and the time interval between 
each train of stimulation was 5 s. The exocytosis was calculated as the in- 
crease in membrane capacitance at the end of eight 500-ms depolarization, 
and the endocytosis was calculated as the decrease in membrane capaci- 
tance at 4 s after depolarization as compared to that at the end of depolariza- 
tion. The decays in membrane capacitance were fitted with either single or 
double exponential functions to extract time constants of the fast and the 
slow endocytosis, as shown previously (He et al., 2008). In Figure SSE, the 
pipette solution contained the following (in mM): 125 Cs-glutamate, 2 Mg- 
ATP, 0.3 Nag-GTP, 1 MgClo, and 0.1 EGTA (pH 7.2 adjusted with CsOH). 
Capacitance increases were evoked by eight trains of depolarization 
(from —70 to 0 mV for 500 ms). 


Optical Tweezer Experiments 

The JPK NanoTracker optical tweezers system was used to measure the 
membrane tether force, which indicates the cell membrane tension. Data 
were collected and calibrated using the PFM Software Release 2.2 equipped 
with the JPK system. To form a tether, a 4-um polystyrene bead approached 
the cell at 1 m/s, arrived at the cell surface, stopped for 10 s for attachment, 
and then was retracted at 1 m/s. 


Imaging Data Analysis 

Alignment of two-color images, semi-automated vesicle-fusion detection, and 
spatiotemporal averaging of local fluorescence signals from multiple events 
were done using a self-written program in MATLAB (MathWorks) (Yuan 
et al., 2015). A 5-pixel (835-nm) diameter circle that encloses the center of a 
fusion spot was selected to calculate the average fluorescence signal of 
VAMP2-pHluorin. LCa-DsRed or Dyni-mKO signals were calculated from 
the same spot, accordingly, and fluorescence intensity changes above 10% 
of the original base value were considered positive for recruitment of the 
marker. The distance between a local maximum VAMP2-pHluorin signal and 
the neighboring LCa punctum was calculated for each fusion event that re- 
cruited endocytic markers. The departure time (t) was defined as the time dif- 
ference between the opening of the fusion pore (as indicated by an increase in 
the fluorescence intensity of pHluorin) and the decrease of the LCa/Dyn1/Syt 
VII puncta signal (Figure 1B). The histograms of t were fitted with a Gaussian 
function (fast population) plus one or two Rayleigh functions (slow populations) 
(Loerke et al., 2009). 

To obtain spatiotemporal averages, all exocytic events were aligned tempo- 
rally according to the frame that exhibited maximum pHluorin fluorescence 
intensity and spatially according to the weighted centroid of the punctum. 
Correspondingly, local endocytic signals (LCa/Dyn1) and the DAG signal 
(PKD-C1ab-GFP) from the other channel were spatiotemporally aligned ac- 
cording to the same time and position. For comparison, the same number of 
signals from the other channel at random frames and sites that was not co- 
localized with exocytosis also were aligned and averaged. 
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To detect CCPs, we first de-noised images with a wavelet transform func- 
tion (Olivo-Marin, 2002), determined positions of puncta by finding local max- 
ima with a mask of 5 x 5 pixels, and then calculated weighted centroids of 
these puncta to yield sub-pixel accuracy. Single-particle tracking was 
accomplished using spatially and temporally global particle assignment, as 
described previously (Jaqaman et al., 2008). Tracks with lengths above ten 
frames were kept and classified as either linear or confined CCP diffusion 
based on a statistical testing algorithm, as previously described (Jagaman 
et al., 2011). Once the trajectory of a CCP intersected within 5 pixels from 
any fusion site, it was considered as possibly related to exocytosis, while 
others were regarded as unrelated to exocytosis. Based on these criteria, 
we compared average velocities of CCPs that were related and unrelated 
to exocytosis (Figure 4H). 

Analysis of microtubule images was performed with Fiji (ImageJ). They were 
first de-convolved using the Richardson-Lucy algorithm (Vonesch and Unser, 
2008), background subtracted (rolling ball radius = 20), and then filtered using 
the steerable plugin (Aguet et al., 2005). 


Statistical Analysis 

All data were analyzed and final images were prepared using Igor Pro software 
(Wavemetrics). Average results are presented as the mean value + SEM from 
the number of experiments indicated. Statistical significance was evaluated 
using either Student’s t test for single Gaussian-distributed datasets or the 
Mann-Whitney rank-sum test for non-single Gaussian-distributed datasets 
(*, **, and *** denote p values < 0.05, 0.01, and 0.001, respectively). 
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